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We can establish a new picture, the perfect fluid strongly coupled quark gluon plasma 
(sQGP) core and the dissipative hadronic corona, of the space-time evolution of produced 
matter in relativistic heavy ion collisions at RHIC. It is also shown that the picture works 
well also in the forward rapidity region through an analysis based on a new class of the 
hydro- kinetic model and is a manifestation of deconfinement. 

1. INTRODUCTION 

Recently, physicists at Brookhaven National Laboratory made an announcement that 
"RHIC serves the perfect liquid" [ H] . The agreement of hydrodynamic predictions [ 
12] of integrated and differential elliptic flow and radial flow patterns with Au+Au data 
at RHIC energies [ El IH El El is one of the main lines of the announcement. We flrst 
study the sensitivity of this conclusion to different hydrodynamic assumptions in the 
hadron phase. It is found that an assumption of chemical equilibrium with neglecting 
viscosity in the hadron phase in hydrodynamic simulations causes accidental reproduction 
of transverse momentum spectra and differential elliptic flow data. From a systematic 
comparison of hydrodynamic results with the experimental data, dissipative effects are 
found to be mandatory in the hadron phase. Therefore, what is discovered at RHIC is 
not only the perfect fluidity of the strongly coupled quark gluon plasma (sQGP) core 
but also its dissipative hadronic corona. Along the lines of these studies, we develop a 
hybrid dynamical model in which a fully three-dimensional hydrodynamic description of 
the QGP phase is followed by a kinetic description of the hadron phase [[7|. We show 
rapidity dependence of elliptic flow from this hybrid model supports the above picture. 
Finally, we argue that this picture is a manifestation of deconflnement transition, namely, 
a rapid increase of entropy density in the vicinity of the QCD critical temperature as 
lattice QCD simulations have been predicted. 



2. SQGP CORE AND DISSIPATIVE HADRONIC CORONA 

A perfect fluid in the QGP phase and the flrst order phase transition to the hadron 
phase is assumed in most hydrodynamic simulations [fSj. While one can flnd various 
assumptions in the hadron phase, e.g. (1) ideal and chemical equilibrium (CE) fluid. 
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(2) ideal and chemically frozen fluid (or partial chemical equilibrium, PCE), or (3) non- 
equilibrium resonance gas via hadronic cascade models (HC). Hydrodynamic results are 
compared with the current differential elliptic flow data, V2{pt), in Fig- 20 in Ref. [|^ with 
putting an emphasis on the difference of assumptions in the hadron phase. The classes 
CE and HC reproduce the pion data well. On the contrary, results from the second 
class, PCE, deviate from these hydrodynamic results and experimental data. In order 
to claim the discovery of perfect fluidity from the agreement of hydrodynamic results 
with V2{pt) data, we need to understand the difference among hydrodynamic results and 
the deviation from data. V2 is roughly proportional to px in low px region for pions. 
In such a case, the slope of V2{pt) can be approximated by V2/{pt)- Integrated V2 is 
generated in the early stage of collisions. Whereas differential f 2 can be sensitive to the 
late hadronic stage since dv2{pT)/dpT ~ ^2/ {pr) indicates interplay between elliptic flow 
and radial flow. In Fig. ^ thermal freezeout temperature T^^ dependence of {px) for 




Figure 1. (pr) for pions as a function of ther- Figure 2. Pseudorapidity dependence 
mal freezeout temperature at 6 = 5 fm at of elliptic flow from hydro and hy- 
RHIC. Solid (dashed) line shows a result with dro+hadronic cascade models are com- 
an assumption of an ideal, chemical equilib- pared with data [Ej. 
rium (chemically frozen) hadronic fluid. 

pions including contribution from resonance decays are calculated at a RHIC energy by 
utilizing hydrodynamic simulations [[HI- {pt) for pions in the chemically frozen hadronic 
fluid decreases with decreasing T*^. This is due to longitudinal pdV work done by fluid 
elements. Whereas {pt) in the chemical equilibrium case increases during expansion. The 
total number of particles in a fluid element decreases due to the assumptions of chemical 
equilibrium and entropy conservation in the hadron phase. Then the total energy is 
distributed among the smaller number of particles as a fluid element expands. These 
are the reasons why the different behavior of {pt) appears according to the assumption 
of chemical equilibrium/freezeout. Under the chemical equilibrium assumption in ideal 
hydrodynamic simulations, increasing [pt) is commonly utilized so far to fix T^^ by fitting 
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Pt slope. However, this is attained only by neglecting data of particle ratio. If particle 
ratios are fixed properly in hydrodynamic simulations to reproduce the data, px slopes, 
especially for protons, are hardly reproduced. The same is true for differential elliptic 
flow: dv2{pT)/PT ~ V2/ {pt) is reproduced by canceling increase behaviors of both t>2 and 
(Pt) under chemical equilibrium assumption [CHI- Agreement of the results from the CE 
model with p^ spectra and V2{pt) data is merely an accident in the sense that this model 
makes full use of neglecting particle ratio to reproduce them. So the HC model turns 
out to be the only model which is able to reproduce particle ratio, pt spectra and V2{pt) 
data. Therefore a picture of the dissipative hadronic corona together with the perfect 
fluid sQGP core is consistent with these experimental data observed at RHIC. 

3. 3D HYDRO AND HADRONIC CASCADE MODEL 

According to the discussion in the previous section, we incorporate a hadronic cascade 
model, JAM [ E], into our previous framework, the "CGC+hydro+jet" model [ IT^. 
Figure 2 shows pseudorapidity dependences of V2 from this hybrid model and ideal 3D 
hydrodynamics with T^^ = 100 and 169 MeV. Here critical temperature and chemical 
freezeout temperature are taken as being = T'^^ = 170 MeV in the hydrodynamic 
model. In the hybrid model, the switching temperature from a hydrodynamic description 
to a kinetic one is taken as Ts„ = 169 MeV. Ideal hydrodynamics with T^^ = 100 MeV 
which is so chosen to generate enough radial flow gives a trapezoidal shape of V2{'r]) [I13j. 
A large deviation between data and the ideal hydrodynamic result is seen especially in 
forward/backward rapidity regions. When hadronic rescattering effects are taken through 
the hadronic cascade model instead of perfect fluid description of the hadron phase, V2 is 
not so generated in the forward region due to the dissipation and, eventually, is consistent 
with the data. So the perfect fluid sQGP core and the dissipative hadronic corona picture 
works well also in the forward region. 

4. SUMMARY: WHAT HAVE WE LEARNED? 

We can establish a new picture of space-time evolution of produced matter from a 
careful comparison of hydrodynamic results with experimental data observed at RHIC. 
What is the physics behind this picture? 7]/s is known to be a good measure to see the 
effect of viscosity where rj is the shear viscosity and s is the entropy density. Figures 3 
and 4 show possible scenarios for temperature dependence of r] and r]/s deduced from 
the discussion in the previous sections. r]/s is small in the QGP phase, which might be 
comparable with the minimum value l/47r [ , and the perfect fluid assumption can 
be valid. While rj/s becomes huge in the hadron phase and the dissipation cannot be 
neglected. Shear viscosities of both phases are found to give ~ 0.1 GeV/fm^ around 
Tc So shear viscosity itself appears to increase with temperature monotonically. 

The "perfect fluid" property of the sQGP is thus not due to a sudden reduction of the 
viscosity at the critical temperature T^, but to a sudden increase of the entropy density 
of QCD matter and is therefore an important signature of deconfinement. 

The author would like to thank M. Gyulassy and Y. Nara for collaboration and fruitful 
discussion. 
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Figure 3. Illustration of the shear viscosity 
as a function of temperature. SYM, HRG, 
and wQGP represent, respectively, super- 
symmetric Yang-Mills model, hadronic res- 
onance gas, and weakly coupled QGP. 



Figure 4. Illustration of the rapid vari- 
ation of the dimensionless ratio of the 
shear viscosity, ri{T), to the entropy den- 
sity, s{T). 
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